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■ Abstract 

Asymptotically flat spacetimes with one Killing vector field are studied. The Killing 
equations are solved asymptotically using polyhomogeneous expansions (i.e. series in pow- 
ers of 1/r and lnr), and solved order by order. The solution to the leading terms of these 
£f~^ • expansions yield the the asymptotic form of the Killing vector field. The possible classes of 

' Killing fields are discussed by analysing their orbits on null infinity. The integrability condi- 

tions of the Killing equations are used to obtain constraints on the components of the Weyl 
CO ' tensor (^o, ^1, VP 2) and on the shear (a). The behaviour of the solutions to the constraint 

equations is studied. It is shown that for Killing fields that are non-supertranslational the 
characteristics of the constraint equations are the orbits of the restriction of the Killing field 
to null infinity. As an application, the particular case of boost-rotation symmetric spacetimes 
is considered. The constraints on ^0 are used to study the behaviour of the coefficients that 
give rise to the Newman-Penrose constants, if the spacetime is non-polyhomogeneous, or the 
logarithmic Newman-Penrose constants, if the spacetime is polyhomogeneous. 



1 Introduction 



The Newman-Penrose (NP) constants |2C(| are a set of five complex quantities, defined for asymp- 



totically flat spacetimes with a smooth null infinity, with the remarkable property of being abso- 
lutely conserved even in the presence of gravitational radiation. Recently ^4j, p6[ , it has been 
shown that in a more general setting — that of polyhomogeneous spacetimes — the NP constants 
are not conserved, but nevertheless, an adequate generalization of them (logarithmic Newman- 
Penrose constants) can be constructed which are indeed conserved. Polyhomogeneous spacetimes 
are spacetimes which are expanded asymptotically in terms of combinations of powers of 1 /r and 
lnr. The introduction of this more general kind of expansion carries a drawback: null infinity 
( J') is no longer smooth. For more details on these, and other aspects of polyhomogeneity, we 
refer the reader to references @, §|, H> and & 

The physical meaning of the NP and logarithmic NP constants is still an open question. 
The work of interpretation has not been easy for a number of reasons. One of them is a lack 
of examples of exact solutions to the field equations representing physically sensible radiative 
asymptotically flat spacetimes. No explicit radiative solution which satisfies all the requirements 
of asymptotic flatness is known. The examples one can make use of reduce essentially to the 
family of boost-rotation symmetric spacetimes [Q , , Q . These spacetimes have two commuting 
Killing vectors, one of them axial, and the other one, such that it leaves invariant the origin's 
null cone. This family includes, among others, the Bonnor-Swaminarayan |Xl") and the C metric. 
Bicak, Hoenselaers & Schmidt have given a systematic method of constructing these spacetimes 

I, i- 
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Boost rotation-symmetric spacetimes describe "uniformly accelerated particles" that approach 
the speed of light asymptotically. The smoothness of the solution requires the spacetimes to be 
reflection symmetric [jjj; therefore at least two particles with opposite acceleration are present, 
and future null infinity contains at least two singular points. The null infinity of a boost-rotation 
symmetric spacetime can be global, in the sense that it admits spherical cuts, but the generators 
are not complete. Ashtekar & Dray || have shown that the C-metric admits a conformal com- 
pletion such that the cuts of JF are the 2-sphere S 2 . This example in particular settled the issue 
of the existence of radiative asymptotically flat solutions to the Einstein-Maxwell field equations. 

The boost-rotation symmetric spacetimes are usually given in a coordinate system that clearly 
exhibits their symmetries. The transformation between this coordinate system and the Bondi 
coordinates used in the asymptotic expansions of the gravitational field has to be given in terms 
of series. These expansions are extremely messy, and usually only the leading terms can be 
calculated explicitly. To add to the problem, the coefficients in terms of which the NP and 
logarithmic NP constants are defined are found deep into the series expansions. As a curious 
observation, Kinnersley & Walker |l7j in a note added in proof mention that the C metric is a 
counterexample to the claim that all algebraically special spacetimes have zero NP constants. 
However, no expression for the conserved quantities is given. 

Bicak & Schmidt Q, starting with a vacuum axially symmetric spacetime, imposed an extra 
symmetry on the spacetime. Expanding the Killing equations in powers of 1/r, and solving 
order by order one finds an asymptotic expression for the new Killing field. The boost-rotation 
symmetry appears playing a privileged role as the only other symmetry an asymptotically flat 
axially symmetric radiative spacetime can have. A generalization of this work for the case of 
electro- vacuum spacetime has been done recently by Bicak & Pravdova poj . 

It is clear that if one wants to pose an initial value problem for a spacetime that has a particular 
symmetry, the initial data cannot be arbitrary, for it should satisfy some constraints imposed by 
the Killing vector field. For instance, if the spacetime is bound to be axisymmetric, the initial 
data cannot depend on (p. In the analysis by Bicak & Schmidt, these constraint equations begin 
to appear when the expansions are carried to higher orders in 1/r. In reference a constraint 
for the news function was found, and in []lo| another for the mass aspect. 

In this article their approach will be generalized in two ways: first it will not be assumed from 
the beginning that the spacetime is axially symmetric, and second it will be assumed that the 
spacetimes can be expanded asymptotically using polyhomogeneous series. The Killing equations 
will be solved using these expansions, and the constraints for the different quantities of interest 
will be deduced from the integrability conditions of the Killing equations. In particular, we will 
be most interested in the constraint on the coefficient of order 6 in 1/r of ^fo, as it is the one 
from which the Newman-Penrose (and the logarithmic NP) constants are calculated. 

As it will be shown, the structure of the constraint equations is very similar to that of the 
continuity equations of the classical mechanics of continuos media. This fact will allow us to gain 
some insight on the behaviour of the solutions to the constraint equations, and put forward a 
tentative interpretation of the physical mechanisms involved, although it may not be possible to 
calculate the explicit form of the solutions. 

The article is organized as follows: in section 2, some preliminaries are discussed. These 
include a brief note on the coordinate system and the null tetrad to be used, the Killing equations 
and their integrability conditions in the NP formalism, some comments on the hypothesis of 
polyhomogeneity, and some remarks on the characteristic initial value that will prove of use on 
later discussions. In section 3 the Killing equations are solved to the first order yielding the 
asymptotic form of the Killing fields compatible with asymptotic flatness. The case of axial 
symmetry is analysed. The result by Bicak & Schmidt Q on the privileged role of the boost- 
rotation symmetry in axisymmetric spacetimes is recovered. Section 4 is devoted to the study 
of the constraint equations that can be deduced from the integrability conditions of the Killing 
equations. The general form of these constraints is discussed, and some general remarks on 
the behaviour of their solutions are put forward. Some emphasis is put on the characteristic 
curves of the differential equations. The constraints for the leading term of the shear (02,0), the 
news function (02,0)1 the mass aspect (Re'l'j' ) an d the coefficients ^q, \I/q, vPq are calculated and 
analysed. The resulting expressions are applied to the particular case of boost-rotation symmetric 
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spacetimes. The final product of this analysis is the constraint equation for the coefficient \& ' 
that gives rise to the NP constants. From these constraints, some remarks on the interpretation 
of these quantities are made. 

There are 4 appendices. Appendices A and B discuss the solutions to two partial differential 
equations that will be used in the body of the article. In appendix C, the orbits of the restriction 
of the Killing field to null infinity are studied, and a classification of them is done. Finally in 
appendix D some spin- weighted spherical harmonics are listed for quick reference. 



2 Preliminaries 

2.1 Coordinates & tetrad 

Most of the calculations will be done with the NP formalism. The coordinates and tetrad used are 
the same as the ones described in Stewart's book |23]. The coordinates x 1 , x 2 , x 3 ) = (u, r, 9, tp) 
are such that u is a retarded time labelling the foliation of the spacetime M. by null hypersurfaces. 
The affine parameter r parametrises the geodesic generators of the null hypersurfaces, and the 
angular coordinates (9, ip) are such that they are constants along the generators of J? and along 
the geodesic generators of the null hypersurfaces. The freedom left in this coordinate construction 
is: a relabelling of the null hypersurfaces, a different choice of angular coordinates, and a rescaling 
and setting of the origin for the affine parameter r. 

The tetrad is constructed so that 1 11 is tangent to the geodesic generators of the null hyper- 
surfaces; is future pointing, and orthogonal to the 2-surfaces u = const. , r = const. (S u ^ r ); 
and m M , and m M span the tangent space to S u ,r (T(5 u . r )). The freedom left in this tetrad con- 
struction consists in: a boost i— > Al^, n^ >— > A -1 ??,^ which yields a rescaling of r; and a spin 
to m i — ► e^mf 1 . 

It can be shown that 

l" = (1) 
= 6$ + Q6£ + C i 8$ 1 (2) 
= e<5f. (3) 

This tetrad is such that e = « = 0, t = W = a + f3, and p and p are real functions. 

The freedom left in the construction of the vectors and gives rise to the notion of 
spin-weighted quantities. The derivatives 8 and 3 work as raising and lowering for the spin weight 
of the relevant quantities. We will stick to Penrose & Rindler's convention |21j]. The connection 
between the 8 and S operators and the directional derivatives 8 and 8 is given by 

3?7 = 8n + s(a — (3)r), (4) 
Bn = 8r) — s(a — (3)<n, (5) 

where 77 is a quantity of spin weight s. A spin- weighted quantity on the sphere (S 2 ) can be ex- 
panded in terms of spin-weighted spherical harmonics. These spherical harmonics are a complete 
and orthonormal set of cigcnfunctions of the operator 39. 

When dealing with differential equations that contain the operators 3 and 3, the following 
lemma will be most useful (2^| . 

Lemma 1 Suppose rj is continuous on the sphere and has spin weight s > 0. Then if 3r/ = 
it follows that rj = 0, and if 877 = then rj is a linear combination of the s Y s ^ m . The analogous 
results hold for s < by interchanging 3 and 3. 

2.2 Poly homogeneity 

In order to handle with ease polyhomogeneous expansions, some conventions will be adopted. We 
will refer to the Newman-Penrose field equations, the Bianchi identities and the frame equations 
in the way they are labeled in Stewart's book ((a) to (r), (Ba) to (Bk), and (Fa) to (Ff)). The 
Killing equations, ( |l6| ) to ( p2| ) will be referred as (Kl) to (K7). A first subscript in the label of a 



3 



given equation will mean that we are just interested in that particular power of 1/r, and a second 
subindex will refer to a particular power of z — lnr. 

Otherwise stated, it will be assumed that all quantities (components of the Weyl tensor, spin 
coefficients, tetrad functions, and coefficients of the Killing fields) are polyhomogeneous functions. 
We will say that a function is polyhomogeneous if in a neighborhood of it can be written as 

oo oo Nk 

/ = £i^" fc = EE^ r ~ felnir > (6) 

k=l k=l j=0 

where fk is a polynomial of degree Nk in z = lnr, whereas fkj is a function of (w, 9, (j)) — no 
dependence on r left — . Let # denote the degree of a polynomial. In the example given above, 
= Nk- Sometimes, the degree of the polynomial will appear next to it in square brackets, 
fk[N k ]. 

The S operator has a polyhomogeneous expansion 

g = r- 1 o i i[0]+r- 2 3 2 [7V 3 + l] + ..., (7) 

where 

B lV (6, tp) = 5^(6, ip) + 2sa x r,{6, p), (8) 

and 5J{r) = 0. 



2.3 The asymptotic characteristic initial value problem 

The asymptotic characteristic initial value problem is usually set by supplying 'J'g on an initial null 
hypersurface Af , a 2 .o on J+ (or and **' , ff on Z = J+C\Nq (or Z Q = J~ n7V ). 

Kannar jl6| has proven the well posedness and existence/uniqueness of the initial value problem 
for data that is C°°. A similar theorem for polyhomogeneous initial data is not yet available. 
However, one can use Kannar's result as a sensible guide when formally solving the initial value 
problem with polyhomogeneous data. 

The most general form for the component of the Weyl tensor 0| is given by 

* = O(r- 3 ln JV3 r), (9) 

that is 

*o = *o[^ 3 ]^ 3 + l-oI^K 4 + $ 5 [N 5 ]r- 5 + ^[iVe]/-" 6 + • • • . (10) 
However, in order to ease the calculations, we will restrict our attention to spacetimes such that 

* = O(r- 4 ln JV4 r), (11) 

This family of polyhomogeneous spacetimes are such that the leading term of the shear (172) 
contains no logarithmic terms (finite shear at J?). Among this family, one finds the "mini- 
mal" polyhomogeneous spacetimes, which are those whose logarithmic terms are directly due to 
non-compliance with the Outgoing Radiation Condition. One has reasons to suspect that these 
particular group of polyhomogeneous spacetimes are the ones with physical relevance p^|,|^|. 
The component 'J'o for these "minimal" spacetimes has the form 

* = + (^'° + mr ) r -5 + ^6,0 + ^6,1 lnr j f -6 + .. . j (12) 

while the non-polyhomogeneous spacetimes obey the Peeling theorem: 

*o = *o'V 5 + *o V 6 + ' ' ' • ( 13 ) 
Details on how to solve the NP hierarchy using polyhomogeneous expansions can be found in 
references pif and |p6| . The results given there will be used whenever it is necessary. 
The logarithmic Newman-Penrose constants are given in terms of "Po by 

Qk = [ *% X hY M )dS, (14) 
Js 2 

where X = maxjiVs, 3A3 + 3, N3 + N4 + 2} p^j . In the case of non-polyhomogeneous spacetimes 
(N3 — N4 — —00, and N5 — so that X = 0) one recovers the original NP constants. 
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2.4 The Killing equations and their integrability conditions in the NP 
formalism 

The Killing equations and their integrability conditions (C^C vp ,\ x = 0, where C vll \ x is the Weyl 
tensor) written in terms of GHP quantities can be found for example in the article by Kolassis 
& Ludwig |Q (a NP version had been worked out previously by Collison & French [Q, but it 
is known to have several typographical errors). The equations given in the present article were 
deduced from theirs, and particularized to the specific NP null tetrad that was described in a 
previous section. 

Let £^ be a Killing vector. It can be written in terms of the vectors of the null tetrad as 

^ = a/" + bn» - cm) 1 - cm?, (15) 

where a, b e K have spin- weight 0, and c e C has spin- weight — 1. The Killing equations (equations 
(Kl) to (K7)) for the tetrad used in this article are: 

Db = 0, (16) 

Aa + (7 + 7)0 = Vc + vc, (17) 

Be = -aa + Xb, (18) 

A6- (7 + 7)6 + Da = 0, (19) 

Sc + cfc = -2pa + 2fib, (20) 

Dc + db = irb-pc-ac, (21) 

Ac -(7 -7)0 + Ba = -To + Ac + Vb + \£c, (22) 

and the integrability conditions for a, ^0, *i and ^3 (equations (IS), (10), (II), (12)) are: 

aDa + bAa - c3er - ccSct + baj - 3ba 7 = 5Q - fQ + (V + 2i,S)o r , (23) 

0.0*0 + oA* - cS^-o - c3* - 4&7*o = 2(7> - iS)$> - 4Q*i, (24) 

aZ)*! + 6A*i - c9*i - c9*i - 267*1 = (V - - Q'^o - 3Q* 2 , (25) 

aD* 2 + 6A* 2 - c9* 2 - c3* 2 = -2Q'*i - 2Q* 3 , (26) 

where 

Q = 36 + 76, (27) 

Q' = 3a - Ifa - Vb, (28) 

V = -Da + Wc + irc, (29) 

2iS = B"c - 3c. (30) 

The equations (|l6|)-(p6|) will be expanded in polyhomogeneous series, and then solved order 
by order. This analysis will yield Killing vector fields that are consistent with the polyhomoge- 
neous asymptotically flat spacetimes of |l3), p4| , p6[ . And more importantly, it will also yield 
constraints for the quantities that arc initial data at No, and Z (i.e. \&0; an d 02,0, *i'° ; ^2'°)- 

It will be assumed that the coefficients in the Killing field can be expanded as polyhomogeneous 
series. 

a = a_ir + a + air -1 + a 2 r -2 + a 3 r -3 + . . . , (31) 
b = b-ir + 6 + bir- 1 + & 2 r~ 2 + b 3 r~ 3 + . . . , (32) 
c = c_ir + c + cir -1 + c 2 r~ 2 + c 3 r~ 3 + . . . , (33) 

where the <Zj, 6^, and are polynomials in z = lnr. 
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3 The asymptotic Killing vector. 



Solving the Killing equations to the leading order in 1 /r yields the asymptotic form of the Killing 
vector fields that are compatible with the asymptotically flat spacetimes under consideration. 
The expansions yield 



d r b 


= o, 


d u a-i 


= o, 


8ic_i 


- o, 


d u b + a_i + a'_ 1 


- o, 


Bic-i + 3ic_i 


= 2a_ 


d-i 


= o, 


d u c-i 


= 0. 



(34) 
(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



From equations (M) an< 4 (PM) together with (k57j) one sees that 



#a_r = #c_! = 0, (41) 

and b is independent of r. Hence, there are no logarithmic terms at this order. From 9ic_i = 
(using lemma 1) one finds that c_i has to be a linear combination of I = 1 spherical harmonics: 

l 

c-i - J2 (- 1 ) m+1 ^- (-iYi,-m) , (42) 

m— — 1 
1 

c-i = 51 3 md y i,m), (43) 

m— — 1 

where A m e C. Using the properties of spin weighted spherical harmonics one can readily find that 

i 1 

0-1 = 2 E {4 + (-l)^- m }(oM. (44) 

m— — 1 

and hence 

l 

6 = ~ J {A m + (-l) m A^ m }( Y hm ) + a(e,ip), (45) 

m— — 1 

where a(0, y?) is an arbitrary integration function. The asymptotic form of the Killing vector field 
is therefore 

/ 1 % esc \ 

C= \-a- 1 u + a(6,<p),a- 1 r,--j=(c- 1 +c- 1 ),- (c_i - c_i) ) . (46) 

Hence, all Killing vector fields compatible with asymptotic flatness can be constructed by provid- 
ing the 3 complex numbers A—i, Aq, A\ and the function a(9, tp). Polyhomogeneity adds nothing 
new at this order. In the axisymmetric case, it has been shown that the integration function 
a{0, if) can be removed with a supertranslation. This suggests that a similar thing can be done 
in the general case. If one lets v! = u + (3, then the Killing vector field will transform as 



e 



(47) 
(48) 
(49) 
(50) 
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In order to remove a, one needs to find a solution (3 to the partial differential equation 

a = c-idifi + c-JS^. (51) 

As is shown in the appendix A, it is always possible to construct a /3 that satisfies equation (fill). 
This shows that the function a(6, if) is associated with supertranslational Killing vector fields, 
i.e. fields of the form 

g ap = (a(6,<p), 0,0,0). (52) 

3.1 Axial symmetry 

In an axially symmetric spacetime, one expects all the quantities to be independent of the co- 
ordinate cp. This fact constrains the expansions in spherical harmonics of functions over S 2 , as 
the only harmonics that are independent of <p are those with m = 0. Hence, in an axisymmetric 
spacetime a quantity of spin weight s, r\ will be of the form 

V = ^2vk(sYk,o). (53) 

k—s 

The Killing vector for axial symmetry is 

rf = (0,0,0,1), (54) 

so that 
and 

a = b = 0. (56) 

Now, we note that the vector can be chosen so that at the origin, its spacelike projection 
lies in a plane containing the axis of symmetry. Due to the axial symmetry of the spacetime, the 
projection will remain in this plane. Therefore n M will have no components in the direction of 
d v . Hence C 3 = 0. The null vectors m M and m M were constructed so that they span T(S u , r )- 
The dimension of T(S Ut r) is 2, but the vector depends on 4 real functions of (u, r, 9). The 2 
extra functions are related to the freedom of performing a spin m* 1 i— > em/ 1 . From m* 1 and 
one can construct two real vectors 

e£ = -J= (m» + m») = (0, 0, Re£ 2 , Re£ 3 ), (57) 

e% = -^( m ^-M») = (0,0,Im^,Im^ 3 ). (58) 

The using the freedom left in the spin it is possible to set Re£ 3 = Im£ 2 = 0. In this way will 
lie in the direction of d v , and £ 2 will be real and £ 3 will be pure imaginary. After this choice, 
all the tetrad freedom will have been removed and all the spin weighted quantities will have lost 
their spin weight. If £ 2 is real, and £ 3 is pure imaginary then from the frame equation (Fb) 

De = P Z? + aT, (59) 

one deduces that a has to be real, and hence from the NP field equation (b) 

Da = 2/W + * , (60) 

is real. These results are summarized in the following lemma. 

Lemma 2 For an axially symmetric spacetime, it is possible to choose the tetrad so that C 3 = 0, 
£ 2 7 a, are real, and £ 3 is pure imaginary. Under this choice, there will be no spin freedom left. 
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3.2 Boost-rotation symmetric spacetimes 



If the spacetime is axisymmetric and happens to possess another Killing vector then 

the two Killing vector fields will form an Abelian algebra . 

If one assumes that the spacetime is axisymmetric, then the coefficients <z_i, b, c_i and c_x 
for the other Killing vector, £ M are constrained to be 

c_i = A sm9, (61) 

c_i =A sin 0, (62) 

a_i = -=(>4o + A )cos6>,' (63) 

b= — p(^ o +^o)cos0. (64) 

Therefore, the asymptotic form of the extra Killing vector field has to be 

^ = (_R e (Ao)u cos 6 + a(8),Re(A )r cos 9,-Re(A) sin 0,-Im(A o )). (65) 

With any loss of generality one can set Im(Ao) equal to zero as the vector (0, 0, 0, lm(Ao)) is just 
a multiple of the axial Killing field. Hence, Bicak & Schmidt's || result has been recovered. The 
resulting vector is known as the boost-rotation Killing vector. 



4 Constraints due to the presence of a Killing field 

4.1 General remarks 

As has been seen in the previous section, the expansion of the Killing equations to order — 1 
in 1/r yields the asymptotic form of the Killing vector field. If one carries the expansions to 
further orders, one will expect to observe the interaction of the symmetry of the spacetime with 
the quantities that are initial data at a given null hypersurface Ao through the integrability 
conditions of the Killing equations. The presence of the Killing vector field will require the initial 
data to satisfy some constraint equations that will be derived from the integrability conditions. 
The constraint equations for 02, and ^| can be deduced from the integrability conditions 

((IS), (10), (II), (12)), equations @, (§f), ©. Looking at the integrability conditions it 

is not difficult to see that the generic form of the constraint equations will be: 

bd u X - c-i^X - c-xdiX + HX = Q, (66) 

where H and F depend on quantities that have been calculated at previous orders in the expan- 
sions. A number of observations can be made from this equation. The constraint equation ( |66|) 
is a linear partial differential equation for the (in principle) complex quantity X. The domain of 
the solutions to this equation is at least a piece of JF =KxS' 2 . This fact suggests that one could 
use separation of variables in order to try to solve the equation when the non-homogeneous term 
Q is not present. Hence, let us assume that 

X(u,9,cp) = U(u)n(6,(p). (67) 
Substitution into the homogeneous part of equation ( |66| ) yields 

(-ua-i)U'Q - c-xUdiQ. - c-xUdiQ, + HUQ = 0. (68) 

Dividing by UCl, and collecting terms one gets 

u> 1 f c-Aii + LAii \ 
~ u 1j-—A n H l m 



s 



Using the classical argument of separation of variables, one sees that the left-hand side depends 
only on u, while the right-hand side depends only on the angular coordinates (9, </>). Hence each 
side has to be equal to a constant A = Ai + 1A2 (Ai, A2eR). Therefore 

U' + AuU = 0, (70) 

and 

c_i5ift + c_icSift- (H + a_ 1 A)fl = (71) 
Equation ( |70| ) can be solved readily yielding 

U{u)=C AlM e-^ + ^>\ (72) 

where C Alt \ 3 is an integration constant. The structure of the function U is extremely suggestive. 
The number Ai is clearly a damping parameter, and A2 is a frequency (which we will assume 
to be positive). In principle there is no mathematical restriction to the values A can take, but 
on physical grounds one would like to restrict them to the first quadrant of the complex plane 
(Ai, A2 > 0) as one expects the "background" gravitational field to damp the propagation of 
gravitational radiation. 

The solution of equation ( |7l| ) will be (in principle) an infinite series of spin- weight s spherical 
harmonics. A general solution for the constraint equation ( |S6| ) is given by integrating over the 
values of the parameter A: 



00 i>00 



X= / C AuA2 e-^+ iA ^ u2 il AuA2 (e,ip)dA 2 dA 1 . (73) 
Jo Jo 

Another interesting observation on the behaviour of solutions to equation (|6^) can be obtained 
from writing the <3i and 3i operators as ordinary partial derivatives. If X is a quantity of spin- 
weight s then 

sin s 9 

(5iX = — — (do — % csc 9d<p) (sin~ s 9X), (74) 
v2 

HiX = Sm 6 (d e + icsc9d v ) (sin s OX), (75) 
v2 

and the differential equation (|66| ) takes the form 

(-a_iu)d u X - ^(c-i + c- 1 )d e X- % °^- 9 (c_ x - c^ 1 )d v X + \ H - - (c_i - c_i)^ X = Q. 

(76) 

Therefore the tangent vector field to the characteristics of equation (U&j is given by 



„ ( 1 t- \ icsc9, \ 
= \ -a-\u + a, -— |(c_i + c_i), -j=-(c„i - c_i) I , (77) 

that is, the restriction of the Killing vector field £ M to J^, Hence, the characteristics 

of equation ( |66| ) are the orbits of The structure of the equations ( |66|) and ( |76| ) is very 

suggestive as it resembles that of the continuity equations of mechanics of continuous media: 
pt + div J — F, where F describes a source or a sink of the quantity p, and J its flux. 

From the theory of first order linear partial differential equations, one knows that equation 
(|76j ) is equivalent to a system of 4 (non-linear) ordinary differential equations: 

dw 

— = -a-iu, (78) 

ft = ->-' +c ->- < 79 > 

icsc9 

d? = -- ^-(^-^o, (80) 

dX / s cot 6* , . \ „ . . 

~dt = -{H-^^-C-^X + Q, (81) 
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with initial data u(0) = u , 9 = f, ip = rj (£e[0,ir], ?ye[0,27r)), X(0) = f(€,rj). The first 3 
equations yield the characteristics/orbits of the restriction of the Killing vector field to . The 
function / is therefore by construction an invariant along the orbits. Note that because we are 
solving the Killing equations and their integrability conditions order by order, then the constraint 
equations will form a hierarchy. It will be necessary to solve all the constraints up to order say 
k in 1/r in order to be able to solve the constraints to order k+ 1. This procedure will give rise 
to a set of functionally-independent invariants along the orbits from which it will be possible to 
construct all the spin-weighted quantities. 

The characteristic curves are represented by the mapping (t, £, rj) i— > (u, 9, tp) (t is the param- 
eter of the curve, and £ and rj are the coordinates of the end point of the curve at the initial cut 
of J") which can be inverted as long as its Jacobian is different from zero. If this is indeed the 



case, one can solve equation (81) to get the solution of the original partial differential equation. 



The solution to the homogeneous part of equation(pl|) can be formally written as: 

X h = /(^) e -/ C ("-^ sc ° te ^- c ->) , (82) 

where £ and rj are functions of (u, 9, tp) as discussed before, and the integration has to be under- 
stood as a line integral along the (unique) characteristic that goes through (u, 9, </>) and which 
can be retrodicted up to the initial cut [Zq = J? C\ A/q). A particular solution to the equation 
( |8l| ) is given by 

X p = e - J c (*->«*«Cs-i-=-i)) J Qe fc(»-^«*0(°-i-c-i)) . (83) 
The complete solution will be therefore 

X = X h + X pi (84) 

where the homogeneous term (Xh) is associated with the propagation of the initial data f(9, <f>), 
which will be damped due to the "interaction" of the quantity X with the background (represented 
by the term (H ~ -^s cot 0(c_i — C-i))X). The non-homogeneous part will account for the added 
effects (in time) of the source/sink. 

4.2 Constraint on the news function 

The shear (a) is a quantity of great physical interest as the derivative of its leading term with 
respect to the retarded time u (&2,o), known as the news function, determines whether or not 
the spacetime is radiative. The outgoing radiation field for the asymptotically flat spacetime is 
determined by the leading term of ^4, which itself depends on the news function (^4'° = — cf). 
In order to obtain a constraint equation for the leading term of the shear (a) one has to expand 
the integrability condition (IS) (equation (p3|)) to order 2 in 1/r. One directly finds that: 

O2,0& - C_lSl(72,0 - C_i3i(72,0 - (O-l + SlC-i - SlC_i)cT2,0 = 0. (85) 

Differentiating the latter with respect to u, one obtains a constraint equation for the news function 
(02,0), 

<72,ob - c_i5i(T2,o - c_i5ict 2 ,o - 2(5ic_i)f7 2 ,o = 0. (86) 

Once these constraint equations have been solved, one can proceed to solve the Killing equations 
@-@ at order in 1/r, 

d - \ (a'_ x + a_x) = 0, (87) 

<3iCo + 3 2 c_i = -CT 2 a-i + A16, (88) 
a = 0, (89) 
3ic + cSiCq + 3 2 c_i + 3 2 c_i = 2a - b, (90) 
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Qib = Co - do - cr 2 c_i, (91) 
c + Sia_i = Aic_i, (92) 

where Ai = <72,Oi an( l ^2 = — 02,o<3i- Note that because of = 0(r _4 ln 4 r), then #02 = 0, 
and hence 02 = 02,0- The coefficients cq and ao can be calculated from equations ( |9l| ) and ( |90| ) 
respectively. It can be easily seen that fj=ao = ^=cq = 0. 

4.2.1 Supertranslational Killing field. 

It is known that an axisymmetric spacetime that admits a supertranslational Killing vector is 
necessarily non-radiative (0 , || . Using equation (|8^) it is easy to show that this assertion is 
still valid even if we remove the hypothesis of axial symmetry. For a supertranslational Killing 
vector field we have that 

c_r = 0, (93) 
a_r = 0, (94) 
b = a(0,cp). (95) 

Therefore the constraint equation for 02.0 is 

a 2fi b = 0. (96) 

This clearly shows that if a ^ then necessarily 02,0 = 0. So we obtain the following 
proposition. 

Proposition 1 If an asymptotically fiat spacetime admits a supertranslational Killing vector field 
then it is non-radiative (&2,o = for all retarded times). 

4.2.2 Boost-rotation symmetry. 

From the discussion in section (2.7), the leading terms of the boost-rotation Killing vector are 
given by 

c-i = ,4osin6>, (97) 

a_i = V2A cos6», (98) 

3ic_i = \/2,4 cos6», (99) 

b = -uV2A cosd, (100) 

where AoeK. If one fixes the tetrad in the way prescribed by lemma 2 in section (2.6), then <7 2j o 
is real, and the constraint equations (|85]) and ( [86] ) simplify to 

(u cot 9)d u (o-2, ) + d e (a 2 ,o) + (cot 6)a 2 , = 0, (101) 

(ucote)d u (& 2 ,o) +d e {& 2 ,o) + (2 cot 6»)ct2,o = 0. (102) 

Both partial differential equations can be solved easily using the method of characteristics. The 
last equation is a particular case of the equation u cot 9d u x + dgx + k cot 9x = H(u, 6) . This kind 
of equation will appear many times more, therefore, a study of its solutions is made in appendix 



B. The general solution of (101) is 



and that of equation(h02J) is 



a 2 .o = (csc0)G( — ), (103) 
u 



1 / sin 6 , 

^2,0 = -a-Fi , (104) 



where G and F\ are arbitrary functions of sin0/u, and F\ — —G'. Hence, we have recovered the 
results of ||. 
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4.3 Constraints on Re ° and ^o ,7V4 

The next step in our study is to obtain constraint equations for Re "J/^' (the mass aspect of the 
spacetime), and VPq (the coefficient that give rise to the logarithmic terms in the expansions). 
The constraint equation for Re ty 3 ,' can be deduced either from the constraint equation (|2^) or 
from the expansions of the Killing equations at order 1 in 1/r: 



-v-ic-x - V2C-1 +ai+ Qicl-x + C%d a a-i + (72 + 72)0-1 = °j 
CT 3 a_i + W 2 a + 3iCi + 3 2 c + 83C-1 = X 2 b, 

I1C1 + (5 2 co + + 81C1 + (5 2 co + 03C-1 — 2[i 2 b + 2psa_i = 2ai, 

a 3 c^i + a 2 c - 2ci + c' : + />;c_i + 3 2 6 - 7r 2 6 = 0, 
Qic_i - icf, + r 2 a_i + 9 2 a_i + 3ia - /U 2 c_i 1 



7C0 + Ci 



+C^9 a c_i - (72 - 72) c -i - A 2 c-i - Aic = 0, 



(105) 
(106) 
(107) 
(108) 
(109) 

(110) 



where 



r- = n- 1M r - 2 + a 3 [N 4 }r- 3 + 0(r~ 4 ln^ 5 r), 
^2 , 1^^4+1 , 



3iMr- 2 + 0(r~ 3 ln JV4+i r), 
4 - Mr" 1 + Ofr" 2 ln^ 1 r 



C 2 "9 a = - [(3i<r 2 ,o)$i + (3i(7 2 ,o)(Si] r- 2 + 0(r- 3 In 



-3 1^^4 + 1 , 



and M = Re*, 



72+72 
^3,0 



-Mr" 2 + 0(r- 3 In 



0, 



(111) 
(112) 
(113) 
(114) 
(115) 



W 2 ' ). From (107) one sees that jj=a\ = 0, and from equation (|109| 



#ci = #(73 = .ZV4. Using all these results equation (105) reads 

— 2M(i_i — C-1B1M + di — 3iCT 2j o3ia-i — 3i(72,o3ia_i — c_i3iM = 0. (H6) 



From equation ( 107 ) one obtains 

a\ = Mb — 3i<72 o3i^ — 3if2 o3i6. 



(117) 



Combining the two equations, and recalling that b = — a_i one finally obtains the desired 
constraint equation, 



bM - c_i3iM - c_i3iM + 3bM = (3iSr 2 , 3i6 + 3icr 2:0 3ifr). 



(118) 



Note that this constraint equation has a sink term ((3i<7 2 ,o3i& + 3i(72,o3ib)) that depends on the 
news function (<7 2 ,o), as one may expect from the Bondi mass formula (see e.g. ^3|). Once M 
is determined, the coefficient ax is readily found. Using the integrability condition (10) (equa- 
tion(p3)) at order 4 in 1/r one can easily deduce the constraint equation for \I/q, 



c^xdi^o + 5-i9i*o + 23ic_x^ = -a_i<; 
in particular the coefficient of the highest lnr power (N4) should satisfy: 

c^x* 4 Q ,N4 + c-i$xK Ni + Zdc-x^ 4 = 0. 



(119) 



(120) 



Note that there are no derivatives with respect to the retarded time, as is a constant of motion 
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4.3.1 Boost rotation symmetry. 



It can be shown that in the case of boost-rotation symmetric spacetimes, the constraint equation 
for the mass aspect M reduces to the one found by Bicak & Pravdova jn| : 

u cot 6d u M + d e M + 3 cot 9M = ud u (d e a 2 fi + 2 cot 0a 2 ,o)- (121) 

From appendix B we learn that the homogeneous part of the solution of the previous equation 
will be of the form Mh = u~ 3 F(u~ 1 sin#). A discussion of the behaviour of the solutions to the 
constraint equation, and their relation with the Bondi mass of the spacetime can be found in }Tq[ ] . 

The constraint equation for ^g' (in a "minimal" polyhomogeneous spacetime) reduces to an 
ordinary differential equation (\l/g is a constant of motion p4j), 



*o' + 2cot6i*o' = 0, (122) 



_d_ 

d<9 

whose solution is 

= Ci csc a 0. (123) 

Therefore 'J'q is singular at 6 = 0, 7r. Note that 'J'o is an invariant of the boost-rotation symmetric 
spacetime, as all the freedom in the tetrad has been removed. Hence, there are singularities at 
the "north" and "south" poles. From this analysis one concludes that a polyhomogeneous boost- 
rotation symmetric spacetime with = 0(r~ 4 ln^ 4 r) can only have a "local" J^, i.e. J' is 
not isomorphic to 5 2 xR §, |. If one wishes to have a boost-rotation symmetric spacetime 
with at least a "piece" of J* {J? ~ 5 2 xM) then one must set 'I' 4 , = 0. Note as well that 
Si^o ° = CiSi(csc 2 6) = 0. This fact will simplify future calculations. 

4.4 Constraint on *}'° and j> 

The integrability condition (II) (equation(p5|)) expanded at order 4 in 1/r yields the constraint 
equation for \E' 1 :0 , while the constraint equation for can be deduced from the integrability 
condition (10) (equation (pi|)) at order 5 in 1/r. In the theory of multipole expansions of stationary 
spacetimes these coefficients are closely related to the dipole moment and the quadrupole moment 
respectively Jl5[ , |l9| . For radiative spacetimes, these relations do not hold anymore; however 
*5>i'° is an indispensable ingredient of all the definitions of angular momentum for asymptotically 
flat spacetimes @, ||. 

The constraint equations are respectively 



b^t - c_i3i¥? - <5_i3i#J - (2o_i + 3ic_i)#J = -o_i*f - 33i6*|, (124) 



and 



- c_!3i^ - c_i3i^ - (a_i + 23 x c_i)^ = 
c 9i*g + c 3i*^ + c_i3 2 ^ + c_i3 2 *o - 43i&*f . (125) 



-a_i^' + (U - a )tf 3' + (4ao - 2b + 2t 2 c_i + 2t 2 c_! + 3ic + 3 2 c_i - 3ic - 3 2 c_i)<J< 4 



Now, from @ one knows that tff = S^ 4 so that (A?4 + l)^ 4 { Na+1 = Si*^ 4 ; therefore the 
only new constraint equation one can deduce from (124) is that for \E' 1 ' (the equations for the 
other coefficients are satisfied identically), therefore: 



bV$'° - c_i3i*i'° - c_i3i^° - (2o_i + 3ic_i)*i'° = -a_i3i*o'° - 33i6*|°. (126) 

This constraint is valid both for polyhomogeneous and non-polyhomogeneous spacetimes. It can 
be regarded as describing a process of transformation of dipole moment (angular momentum) 
into mass monopole momment due to the gravitational radiative process. 

For a non-polyhomogeneous spacetime the constraint for the coefficient "Jq' reduces to, 

&*o'° - c-i8i*o'° - c_i3i^° - (o_i + 23ic_i)^'° = -49x6* (127) 
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while for a "minimal" polyhomogeneous spacetime the leading coefficient (^g' ) should satisfy 

b^ 5 A - c-iSi^ 1 - c^d^ 1 - (a_i + 25ic_ 1 )^' 1 = -4Si6Si*J°. (128) 

In a similar way to what happened with the constraint equation for vf^' , these two last equations 
can be interpreted as describing a process of interchange of mass quadrupole moment into dipole 
moment (angular momentum). Note that for the "minimal" polyhomogeneous spacetime there 
will be two coefficients associated with the quadrupole, the logarithmic one being dominant far 
away from the source. 

As in the previous subsections, once the constraint equations have been solved, one can easily 
find the coefficients a 2 and c 2 from the Killing equations (K5) and (K6) (equations (^0|) and (pl|)) 
at order 2 in 1/r: 

S3C0 + 9 2 ci + Sic 2 + 04C-1 + 2p 3 a + 2p 4 a_i - 2/i 3 6 

+34C_i + 9 3 c + 5 2 ci + 9ic 2 = 2a 2 , (129) 

p 3 c + P4C-1 - 3c 2 + 3 3 6 + c' 2 - 7r 3 6 + ct 2i0 ci + a^c-i + CT3C0 = 0. (130) 
From these equations one can see that #a 2 = ^73 = A4 + 1 and #c 2 = A5. 



4.4.1 Boost rotation symmetry. 

The constraint equations for VE"^' are t 
mogeneus boost-rotation symmetric spacetimcs: 



The constraint equations for VE"^' are the same for non-polyhomogeneous and "minimal" polyho- 



u cot 6d u vf° + d e ^t'° + 3 cot = 3utf ij' . (131) 

The constraint equation for the leading coefficients of for a non-polyhomogeneous spacetime 
(^g' ), and a "minimal" polyhomogeneous spacetime ('I'g' 1 ) are respectively, 

u cot 9d u ^ 5 ' + d e ^Q° + 3 cot 6^'° = 2u#t'°, (132) 

and 

u cot edu^l' 1 + de^ 1 + 3 cot 6^ 5 A = 0. (133) 

The homogeneous part of the solutions of all these equations will be of the form u~ 3 F(u~ 1 sin 9). 

Note that in the "minimal" polyhomogeneous case, no source/sink term occurs. The solution 
in this case is therefore: 

ei 1 „ / sin 6 \ , 

*o' = — *2 > 134 

w J \ u J 

where F2 is an arbitrary function of the argument. Now, using the evolution equation for ^q' 
(see or H|), one finds that 

= B^f 1 = 0, (135) 
whence VPg' 1 will also be a constant of motion for boost-rotation symmetric spacetimes. Hence 

to' 1 = C 2 csc 3 <9. (136) 

Again one can see that poly homogeneity in these class of spacetimes gives rise to a local J? . 

4.5 Constraints on \&q and the NP constants. 

Finally, we are able to deduce constraint equations for As seen in section 2.3, the logarithmic 
NP constants are given in terms of an integral of ^S> ' , where X = A5 if A3 = —00. The 
constraint equation will be much more complicated in this case. The expansion at order 6 in 1/r 
of the integrability condition (10) (equation (p4[)) gives an equation of the form: 

6*[j - c-tfitfg - c-i9i^ - 2(o-i + Sic-i)* 6 , = K, (137) 
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where K is a complicated expression depending on ^g, 'J'q and their Si and Si derivatives, and 
on \&i'°, ^2'' °0j a i> 6, & c -1j c 0j c i- In the particular case of the non-polyhomogeneous 
spacetime the coefficient which yields the Newman-Penrose constants ^q' satisfies, 



- c-iSi^g' - c_iSi<° - 2(o_i + 3ic_i)*g' = 
(§6 + 8c-i3ict 2 ,o + 6c_i3icr 2 ,o - 2Sic_icr 2 ,o + 23ic_ict 2 ,o)^'o' 

-c_ia 2 ,oSi*o'° + c 3i*o'° - c_ia 2 , 3i*o'° + So3i*o'° 

-49x6*5'° + 4cr 2 ,o3x6*^ - 4&r 2 *^°, (138) 

and the analogous coefficient (^g' 1 ) for a "minimal" polyhomogeneous spacetime has to satisfy a 
similar constraint equation 

b^' 1 - C-iSi^' 1 - C-xBx^' 1 - 2(o_i + Sxc-x)^' 1 = 
(|6 + 8c_i3ict 2i o + 6c_i3icr 2: o - 2Sic_ict 2i o + 23ic_ict 2i o) v I / o' 1 
-c-ict 2 ,oSi*o' 1 + coSi^o' 1 _ c_ict 2 ,oSi*o 1 + coSi^o' 1 

+43x63x*^ + 4 f 7 2 ,oHx63x*o'° - 46t 2 3i^'°, (139) 

The source/sink terms in these two equations are much more complicated than those in previous 
sections, and hence their interpretation is not that clear cut. 

Using the ideas and notation of section 4.1, one can split into its homogeneous (^n' )h 

and non-homogeneous parts (4 , o' X ) p . From equations (|8^) and ( |S3"| ) one finds that (as path of 
the line integrals reduces to a point) 

(*l' x )hWo = M<p), (wo) 

and that 

(<' X )pWo = 0. (141) 

The Newman-Penrose constants can be evaluated on any null hypersurface, in particular at A/"o, 
therefore 

Qk = [ f(e,v)( 2 Y 2 , k )dS. (142) 
Js 2 

The conservation of Q? shows that the added effect over time of the complicated source terms 



of equations ( |138| ) and (|139| ) cancels exactly the damping of the initial data f(9, < 



4.5.1 Boost rotation symmetry. 

The constraint 
are of the form: 



The constraint equations for ^q' and vE'q' 1 in the case of boost-rotation symmetry spacetimes 



ucotffflu^' ) + d e(^a a ) +4cot0(*[j' ) = H x . (143) 

For the "minimal" polyhomogeneous spacetime the non-homogeneous term simplifies to (recall 
that ct 2 . = u~' 2 F\{u~' 1 sin0) and ^ 5 A = C 2 esc 3 9) 

( sin 9 \ 

Hi = -C 2 csc 3 (6») cot(6»)u- 2 Fi f — J . (144) 
Finally, the homogeneous part of the solutions to the constraint equations have the form 

Aft ■ (»5) 
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5 Conclusions 



We have seen that if an asymptotically flat spacetime — polyhomogeneous or not — is assumed to 
have a Killing vector field, then using the integrability conditions of the Killing equations it is pos- 
sible to obtain constraints on the different components of the Weyl tensor and the news function. 
These extra equations together with the evolution equations derived from the Bianchi identities 
suggest the existence, in non-stationary spacetimes, of processes of transformation of multipole 
moments of a given kind into others of a different class. Unfortunately, this interpretation is done 
in terms of quantities that can only be defined rigorously for stationary spacetimes. 

Newman & Penrose [|o) found that for a stationary spacetime the Newman-Penrose constants 
have the structure 

(dipole) 2 — (monopole) x (quadrupole). (146) 

If one considers a system that initially is stationary, and later undergoes a process of gravitational 
radiation, finally settling down into a stationary state, one can see that the NP constants impose 
a "selection rule" for the class of final states achievable. The idea of transformation of multipole 
moments discussed above fits with this idea, as the monopole moment of the source will be 
radiated according to Bondi's mass loss formula, the dipole and quadrupole moments changing 
accordingly in order to preserve the value of the NP constants. 
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A The solution to the equation: a = c -\d\f3 + c -\<5\fi. 

The function (3 is of spin weight zero, hence 

w -£{!-ib£}* (147 » 

Write c_i = c r + ici. Then we see that 

C _ 1 9 1 /3 + c_ x 9 1( 3 = V2 (cr% + -^^) , (148) 
\ o0 sin dtp J 

hence the partial differential equation to be solved is 

* = ^f + ^#Y (149) 



d0 sin 9 dtp 

Recall that both a and (3 can be expanded in terms of spherical harmonics as 



a = ttlmYlm, (150) 

I m 

= ^^Am« m . (151) 



I m 



Now, the right-hand side of equation ( |149| ) can also be expanded in terms of spin zero spherical 
harmonics as 
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LAI, 



(152) 



L M 



where /lm is linear in the /3j m . This gives rise to the following infinite set of linear equations in 

OiLM = fLM{A m ,0lm), (153) 

that in principle can be solved to any desired order, yielding a solution to the equation (|149| ). 



B The solution to the equation: ucot9d u x + dgx + k cot Ox = 

H{u,e). 

It is not very complicated to find the solution to the initial value problem 



dx dx 

u cot 9— — h + k cot Ox — H(u, 9), 
ou oO 



with initial data 



x{l,6) = F{6). 
The associated system of ordinary differential equations is 



(154) 



(155) 



du 
drj 
cW 
drj 
dx 
di] 



ucotO, 



-kcotx + H(u, 9), 



(156) 
(157) 
(158) 



with initial data 



Equation (157) gives 



u(0) = 1, 
0(0) = t 
i(0) = F(0« 



t + v, 



and so we can solve now equation (|15^) yielding 



sin^ 



sin(?7 + f). 



Therefore the characteristics of the differential equation are given by 



1 



sin£ 



■ sin f 



(159) 
(160) 
(161) 



(162) 



(163) 



(164) 



It will be necessary to invert this last expression in order to have £ as a function of u and 9. The 
range of 9 is in principle [0, it). There will be problems with the invertibility whenever u = sin#. 
The solution will break down there. This comes from the fact that the polar coordinates (9, <p) 
are not good coor dinates for the sphere, and several coordinate patches are needed to cover it. 
Equation (158) is more involved, 
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dx ( 1 \ 

— = -k cot(£ + r/)x + if -— sin(?7 + £)> C + »7 • (165) 

(277 \ Sm S ' / 

The solution to the homogeneous equation can be obtained directly by integration 

1 /sin0\ 

x h = — F \ » 166 

u k \ u J 

while the particular solution to the non-homogeneous equation is given in an integral form 



x p = csc fc (77 + £) 



if 



sin (t + £) 



10 \ sin £ 
Then the solution to the differential equation is 



,t + n sin* (t + £) dt 



(167) 



x(u, 9) 



1 /sinf 

— F ( 

u \ u 



■ csc k {9) 



(168) 



- ^— arcsin(sin ^/u) 



The key to interpreting the solution (|168|) is to regard equation (|154| ) as a continuity equation 
completely analogous to the p, t +V • J = Source continuity equation of fluid mechanics. The 
terms ucotOd u x + dgx correspond to the time derivative plus divergence bit of the continuity 
equation, while the k cot 9x term is a damping term due to the interaction of the gravitational 
field with itself (note that the constant k gives rise in the solution to a 1 /u k term that diminishes 
the amplitude of the initial data. Finally, the function H works like a source/sink term. 



C The orbits on the sphere of the Killing vector fields 

As seen in section 4, the restriction of the Killing vector field to is given by 



-a-iu + a, 



= (c_i+c_i),- 



-(c-i 



c-i. 



(169) 



V2 V " V2 

The integral curves of this vector field in £ can be conveniently visualized as curves on the unit 
sphere parametrized by the retarded time u. The structure of this Killing vector field is not very 
complicated, depending only on 3 spherical harmonics. Therefore, it is quite tempting to attempt 
a study of these orbits. 

Throughout this article, we have been using the angular coordinates (9, (p) on the cuts of 
J^, however the study of the orbits is more easily done using the stereographic coordinates 
= cot ^9e te ), and working on the complex plane plus the point at infinity. The iYi m spherical 
harmonics in stereographic coordinates are given by p3[ : 



■ y " " -ir«TTC ( (170) 
' y '° " -V^TT^ (171) 

Now, the critical points of the vector field correspond to the points on the sphere (complex plane) 
where c_x = 0. Using the spin- weighted spherical harmonics in spherical coordinates one can 
prove the following lemma. 
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Lemma 3 The restriction of a Killing vector field of an asymptotically flat spacetime to J? 
vanishes at most in two points of the sphere (and at least in one). 



The result follows noting that if c_i = 0, then one has 

A^+AoC + A^ = 0, (173) 

a second degree equation in C. From the Fundamental Theorem of Algebra, one knows that the 
equation has two roots. 

The vector field will vanish at least in one point, so without loss of generality one can set this 
point to be the origin of the complex plane (the North Pole of the sphere) , and the other root (if 
present) to lie on the real axis. Hence one can write 

- Ci = e i.^ + K (174) 
i + CC 

Let £ = x + iy, then 



(a(x 2 — y 2 ) + bx) cos a; — (2axy + by) sin a; 
1 + x 2 + y 

(a(x 2 — y 2 ) + bx) sin a; + (2axy + by) cos uj 
1 + x 2 + y 

And the orbits on the complex plane are given by 

dy Imc_i (a(x 2 — y 2 ) + bx) sincj + (2axy + by) cosw 



_ —y T LUSH/ - ililj T Olliu/ 

Rec_i = — — — — - , (175) 

j m ~ _ —y } T UJU) Bill UJ + (2u,juy t uyj ixjs w (176) 



(177) 



da; Rec_i (a(x 2 — y 2 ) + bx) cos uj — (flaxy + by) sinw 

There is no integrating factor for this equation for an arbitrary lj, however the form of the 
orbits can be found readily by noting that they will have horizontal tangent at the "rotated" 
hyperbola 

asinwx 2 + 2acosuixy — asinujy 2 + fesinw + b cos uiy = 0, (178) 
and vertical tangency at another hyperbola 

acoswx 2 — 2as'muixy — acosuiy 2 + bcoscux — bsinujy = 0. (179) 

These two curves intersect by construction only at (0,0) and (— b/a, 0) (the critical points). If 
tt/2, 7r, or 37r/2 then none of the hyperbolae are degenerate (i.e. they are not intersecting 



lines). The degenerate cases will be studied separately. Linearizing the differential equation (177) 
around the origin one obtains (if 6^0) 

dy sinujx + cos ujy 

-7- — : , (J-oUJ 

da; cos ujx — sin ojy 

The solution of this equation is the logarithmic spiral 

x 2 _|_ y2 _ Q e 2cotuarctm(y/x) (181) 

A similar result follows when linearizing around (— b/a, 0). Therefore if a and b are different from 
zero and u ^ 0, 7r/2, tt, or 3tt/2 then the orbits spiral around (0, 0) and (—b/a, 0). 



In the case that uj = 0, tt then the differential equation for the orbits (177) can be readily 
solved yielding 

That is, the orbits are circles that intersect each other in the points (0,0) and (— b/a, 0). 

Now, if a, b ^ and ui — tt/2, 3ir/2 then the orbits are non-intersecting circles whose centres 
lie on the x axis, and have a reflexion symmetry along the line x = —b/2a. 
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Figure 1: Spiral orbits on the complex plane and on the sphere. 




Figure 2: "Electric dipole" -like orbits on the complex plane and on the sphere 



If b = then the two critical points coincide, and one gets a family of tangent circles at the 
origin. The angle uj yields the orientation of the family. 

If a = 0, then one of the critical points is at the origin (North Pole), and the other at infinity 
(South Pole). The generic orbit (uj ^ 0, tt/2, it, 3tt/2) is the logarithmic spiral described by 
equation (181). If u — 0, ir then one obtains lines that intersect at the origin and at infinity, 

y = Cx. (183) 

while if to = 7r/2, 37r/2 then the orbits are concentric circles. 

The previous discussion gives a complete classification of asymptotically flat spacetimes with 
one Killing vector field in terms of their orbits at <f , The results are collected in the following 
proposition: 

Proposition 2 Let M be an asymptotically flat spacetime with complete J? and one non super- 
translational Killing vector field. Then the orbits of the restriction of the Killing field to J? (£f 1 \j ) 
regarded as curves on S 2 parametrised by u satisfy one of the following: 

1. there are two critical points, and the orbits spiral around them; 

2. there are two critical points, and the orbits connect the two points ("electric dipole" struc- 
ture ); 

3. there are two critical points, and the orbits are periodic around them; 

4- there is one critical point, and the orbits are closed curves tangent at them. 

Observe that the critical points are not necessarily antipodes (i.e. located at opposite points on 
the sphere). 



20 



Figure 3: Orbits with a single point on the complex plane and on the sphere. 

y 





Figure 4: Periodic orbits on the complex plane and on the sphere. 

D Some spherical harmonics 

D.0.2 Spin-weight 



D.0.3 Spin-weight 1 



5*0,0 



47T 



0*1,1 = -\/^sin0e^ 

»7T 



0*1,0 
0*1,-1 



_3 

47T 

/]T 

8^ 



cos 9 
sin6»e-^ 



1*1,1 

1*1,0 

i*Vi 



— (cos. + IK 



sm( 



— (cos.- IK 



i(f> 



(184) 
(185) 
(186) 
(187) 



(188) 
(189) 
(190) 
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D.0.4 Spin-weight 2 



2^2.2 


" 3 \ 


/-|-(l + cos0) 2 e 2 ^ 


(191) 


2*2,1 




/— sinflfl + cosflW* 

247T 


(192) 


2*2,0 


CO 1 CQ 


/-^sin 2 

V 47T 


(193) 


2^2,-1 


= 3 \ 


/-^-sin^l-cos^e - '* 

247T 


(194) 


2*2,-2 


= 3 \ 




(195) 
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